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We investigate positive-parity states of with the calculation of antisymmetrized molecular 
dynamics focusing on pn pair correlations. We discuss effects of the spin-orbit interaction on energy 
spectra and pn correlations of the J^T = = 35^0, and states. The state has almost 

no energy gain of the spin-orbit interaction, whereas the 3]*’0 state gains the spin-orbit interaction 
energy largely to come down to the ground state. We interpret a part of the two-body spin-orbit 
interaction in the adopted effective interactions as a contribution of the genuine NNN force, and 
find it to be essential for the level ordering of the 3]''0 and 1]''0 states in ^°B. We also apply a 2a-|-pn 
model to discuss effects of the spin-orbit interaction on T = 0 and T — 1 pn pairs around the 2 q; 
core. In the spin-aligned J^T = 3^0 state, the spin-orbit interaction affects the {ST) = (10) pair 
attractively and keeps the pair close to the core, whereas, in the 1+0 state, it gives a minor effect 
to the (ST) = (10) pair. In the 0+1 state, the (ST) = (01) pair is somewhat dissociated by the 
spin-orbit interaction. 


I. INTRODUCTION 

In the progress of experimental researches on proton- 
rich nuclei, the interest of proton and neutron {pn) pair 
correlations has been revived in these years. In the study 
of pn pairing, the competition between isoscalar T = 0 
pairing and isovector T = 1 pairing is one of the essential 
problems m Z N nuclei The nuclear interaction 

in a free space is more attractive in the T = 0 spin-triplet 
even (^E) channel than in the T = 1 spin-singlet even 
(^E) channel as known from the bound state, deuteron, 
formed by two nucleons in the ^E channel. However, at 
the nuclear surface and in nuclear medium, the competi¬ 
tion between T = 0 and T = 1 pn pairs occurs. Because 
of the stronger ^E interaction than the ^E interaction, it 
is naively expected that the deuteron-like T = 0 pair is 
more favored than the T = 1 pair as seen in the ground 
state spin, J'^T = 1+0, of ®Li and ^®F. However, the 
T = 1 pair is often favored rather than the T = 0 pair in 
medium- and heavy-mass regions as seen in the ground 
state spins of Z = =odd nuclei because the spin-orbit 
mean potential favors the T = 1 pair m- Moreover, the 
spin-orbit potential favors a spin-aligned T = 0 pn pair 
[12 - [l^ . These facts indicate that the spin-orbit interac¬ 
tion plays an important role in the competition between 
T = 0 and T = 1 pn pairs in nuclear systems. 

Investigations of Z = = odd nuclei are helpful to 

understand features of pn pairs at the nuclear surface. 
Based on a three-body picture of a core nucleus with 
two valence nucleons, one can discuss the competition 
between T = 0 and T = 1 pn pairs from the ordering of 
jTTT — 2 ^-i-Q Q-i-j^ states. For example, ®Li and 
have the J^T = 1+0 ground states and the J^T = 0+1 
excited states indicating that the T = 0 pair is favored 
rather than the T = 1 pair. On the other hand, in "'^^Sc, 
the ground state is J'^T = 0+1 because the T = 1 pair is 
favored by the spin-orbit potential at the surface of the 


^°Ca core as discussed by Tanimura et al. based on a 
three-body model calculation [l^. 

In the previous paper [T^, two of the authors, Kanada- 
En’yo and Kobayashi, discussed effects of the spin-orbit 
interaction on pn pairs at the surface of in ^®F based 
on an ^^0-\-pn model, and found that the level structure 
of J'^T = 1+0, 0+1, and 3+0 states is affected by the 
strength of the spin-orbit interaction. Namely, the spin- 
orbit interaction reduces the T = 1 pair energy in the 
0+1 state, and it largely contributes to the energy of a 
spin-aligned T — 0 pn pair attractively to lower the 3+0 
energy, whereas it gives a minor effect to the T = 0 pair 
energy in the 1+0 state. 

In ^°B, the ground state is the 3+0 state and the first 
excited state is the 1+0 state at = 0.72 MeV. Based 
on a 2a-\-pn picture, this fact indicates that ^°B is an 
interesting system in which the level inversion between 
the 1+0 state having a T = 0 pair in an S'-wave (a 
pair moving in the total-angular-momentum L = 0 state 
around the core) and the 3+0 state having a spin-aligned 
T = 0 pair (a pair moving in a total-angular-momentum 
L — 2{D-wave) state around the core) occurs. In these 
years, ab initio calculations using the no-core shell model 
(NCSM) approach with realistic nuclear forces based on 
the chiral perturbation theory [13 have been achieved for 
H ^ 10 nuclei [l^. The NCSM calculations with effec¬ 
tive interactions derived from the chiral two-body {NN) 
and three-body {NNN) nuclear forces [l^, [2^ describe 
well the experimental low-lying spectra of ^*^B and show 
that the NNN force is essential to reproduce the order¬ 
ing of the 3+0 and 1+0 states in ^°B. Recently, Kohno 
pointed out that the NNN force provides an attractive 
contribution to the effective two-body spin-orbit interac¬ 
tion in nuclear medium (2lj| . Therefore, it is expected 
that the NNN force may also contribute to pn pairs in 
nuclei through the effective spin-orbit interaction. 

In this paper, we investigate structure of ^°B and clar¬ 
ify effects of the spin-orbit interaction on T = 0 and 
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T = 1 pn pairs based on the calculation of antisym¬ 
metrized molecular dynamics (AMD) [22l - [^ using phe¬ 
nomenological effective nuclear interactions. The AMD 
method is a model for structure studies and has been 
proved to be one of the successful methods for light nu¬ 
clei, in particular, to describe cluster structures of ground 
and excited states. For instance, 2a cluster structures 
of neutron-rich Be isotopes are described systematically 
with the AMD calculations [l^, [2^. We calculate 
with the AMD method and find that 2a cluster core is 
formed in ^°B. We discuss the role of the spin-orbit inter¬ 
action in energy spectra and features of a pn pair around 
the 2a core in ^°B. Moreover, we associate a part of the 
effective two-body spin-orbit interaction with the NNN 
force based on the Kohno’s evaluation, and discuss its 
effect on the ^°B energy spectra. We also discuss ®Li 
spectra having a pn pair around an a core for compari¬ 
son. 

This paper is organized as follows. In Sec. im we ex¬ 
plain the formulation of the AMD method. In Sec. Illli 
the calculated results for ^°B are shown. We discuss ef¬ 
fects of the spin-orbit interaction on the ^°B energy spec¬ 
tra based on the AMD result in Sec. El In Sec. El we 
perform an analysis using a 2a -I- pn model to discuss ef¬ 
fects of the spin-orbit interaction on the pn pair around 
the 2a core. A summary is given in Sec. IVII 


II. FORMULATION OF AMD AND EFFECTIVE 
NUCLEAR INTERACTIONS 

A. AMD method 

We apply the method of the variation after parity and 
total-angular-momentum projections of the AMD model 
(AMD- 1 -VAP) [ 2^,[23 to obtain A-nucleon wave functions 
for the ground and excited states of a nucleus with the 
mass number A. We here briefly explain the formulation 
of the present AMD calculation. 

An AMD wave function is given by a Slater determi¬ 
nant. 


d>AMD(Z) = 


-^A{v3i,</52,...,(^a}, 


( 1 ) 


where A is the antisymmetrizer, and the ith single¬ 
particle wave function is written by a product of 
spatial((/>i), intrinsic spin(xi) and isospin(ri) wave func¬ 
tions as 


‘P^ = (2) 

/ o,,\ 4/3 -v-. 

0x.(rj) = j exp{-v{rj - (3) 

Xi = (4) 

(()Xi and Xi are spatial and spin functions, respectively, 
and Ti is the isospin function fixed to be up (proton) 


or down (neutron). Accordingly, an AMD wave func¬ 
tion is expressed by a set of variational parameters, 
Z = {Xi,X 2 ,..., Xa,Ci,6, ... ,Ca}, which specify cen¬ 
troids of single-nucleon Gaussian wave packets and spin 
orientations for all nucleons. 

The parameters Z are determined by the energy vari¬ 
ation after parity and total-angular-momentum projec¬ 
tions to obtain the optimized AMD wave function. 
Namely, in the AMD-1-VAP method, X^ and = 1 ~ 
A) for the lowest J'^ state are determined so as to min¬ 
imize the energy expectation value of the Hamiltonian, 
(<I>|i7|<I>)/($|$), for the J'" eigen wave function projected 
from the AMD wave function; $ = P^’^$amd(Z). Here, 
^MK ^^le parity and total-angular-momentum projec¬ 
tion operator. For each J'^ state, the optimum set Zj^ 
of parameters is obtained. After the VAP, to describe 
states, we superpose the J’^-projected AMD wave func¬ 
tions expressed by the obtained parameter sets Zj,^^, for 
various J''" as 

4^4- = E (5) 

J'^' ,K 

where coefficients aj^{J'TT'] K) are determined by di¬ 
agonalizing the norm and Hamiltonian matrices. For 
Z = N = odd nuclei, T = 0 and T — 1 projections 
are approximately done by using the proton-neutron ex¬ 
changing operator Pp^^n after the energy variation as 

= E {ajj(^V;A) + 6jj(jV;iX)Pp^„} 

J'^\K 

4’AMD(Zy^,), (6) 

where aj^ {J'tt'; K) and ( J'tt'; AT) are determined by 
the diagonalization. For T = 0 and T = 1 states, 
aj^{J'TT'-K) Ri -6jj(J'7r'; AT) and ajj (J'tt'; AT) r; 

6jj(J'7r'; AT) are obtained, respectively. In the present 
framework, we do not explicitly assume ojj ( J'tt'; AT) = 
±6,/j(J'7r'; AT) because the isospin symmetry is slightly 
broken in the Hamiltonian because of the Coulomb force. 
However, the obtained 'kjj for states of ^°B and ®Li 
are found to be approximately T eigen states and can be 
assigned to experimental J'^T states. 

For the width parameter v of single-nucleon Gaussian 
wave packets, we choose v = 0.235 fm“^ which is the 
same value used for and ^^B in Ref. [1^ and was 
originally determined for ®Be in Ref. [ 2 ^. 

In the AMD framework, existence of clusters is not 
assumed a priori because Gaussian centroids X^ of all 
single-nucleon wave packets are independently treated as 
variational parameters. Nevertheless, if the system ener¬ 
getically favors a specific cluster structure such the struc¬ 
ture is obtained in the energy variation because the AMD 
model space contains wave functions for various cluster 
structures. Therefore, the AMD method is suitable to 
investigate whether the clusters are formed or not in the 
system. 
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Note that the AMD wave function is similar to the 
wave function used in Fermionic molecular dynamics cal¬ 
culations Bill, though some differences exist in width 
parameters of single-nucleon Gaussian wave packets and 
the variational procedure. Another difference in the 
AMD and FMD calculations is effective nuclear inter¬ 
action. In the AMD calculations, phenomenological ef¬ 
fective interactions are usually used differently from the 
recent FMD calculations, in which effective interactions 
constructed from the realistic nuclear force by means of 
the unitary correlation operator method are used [3l|. 


B. Effective nuclear interactions 

We use the finite-range central and spin-orbit interac¬ 
tions as effective two-body nuclear interactions, 


12 = Vcir){w + bPcr - hP-r - mPcrPr) 
+ I’lAr )—;;- 


(7) 


where Pa- and Pr are the spin and isospin exchange op¬ 
erators, r is the relative distance r = |r| for the relative 
coordinate r = ri — r 2,1 is the angular momentum for r, 
and s is the sum of nucleon spins s = -|- S 2 . We ignore 

the term of the spin-orbit interaction. In the present 
paper, we use the Volkov No.2 central interaction 32|, 


For the strengths of the spin-orbit interaction, we take 
uis = ui = —U 2 - Uls is the strength parameter of the 
effective spin-orbit interaction and, in principle, it may 
depend on nuclear systems reflecting contributions from 
the three-body force and the tensor force as well as the 
original spin-orbit force in bare nuclear forces. It may 
also have structure model dependence, and therefore, is 
considered to be an adjustable parameter in model calcu¬ 
lations. In the present paper, we use uis = 1300 MeV to 
reproduce the Is splitting between 3/2“ and 1/2“ states 
in ^Be in the AMD- 1 -VAP calculation. We also use a 
slightly weaker strength uis = 1000 MeV to see the de¬ 
pendence of energy spectra on the strength uis of the 
spin-orbit interaction. The strength of the effective spin- 
orbit interaction can be estimated by the Scheerbaum 
factor Bs MM defined as 

2 tt r°° 

Bs = -/ drr^ji{qr)vis{r), (10) 

9 Jo 

with q = 0.7 fm“^. Here ji is the spherical Bessel func¬ 
tion. For the G3RS spin-orbit interaction with uis = 
1300 MeV and 1000 MeV, Bs equals to 103 MeV and 79 
MeV. 

In Table HI we list the adopted interaction parameter 
sets of effective nuclear interactions labeled (A) and (B) 
with the strength uis = 1300 MeV and (A’) and (B’) 
with Uls = 1000 MeV. 



with vi = —60.65 MeV, V 2 = 61.14 MeV, oi = 1.80 fm, 
and 02 = 1.01 fm, and the G3RS spin-orbit interaction 


vis{r) = Ml exp 



-I- M 2 exp 



with bi = 0.60 fm and 62 = 0.447 fm. 

For the Volkov central interaction, we use the Wigner 
and Majorana parameters, w = 0.40 and m = 0.60, which 
reproduce the a-a scattering phase shift, and the Bartlett 
and Heisenberg parameters, b — h = 0.125 which repro¬ 
duce the deuteron binding energy. The b and h are the 
parameters which can control the ratio / of the in¬ 
teraction to the ^E interaction for the fixed w + m value 
as / = (w + m + b + h)/{w + m — b — h). The ratio is 
/ = 1.67 for the present parametrization. Generally, in 
effective two-body central interactions for structure mod¬ 
els, the ratio may change depending on nuclear systems 
because of medium effects and it is usually somewhat 
suppressed in nuclei. Therefore, b and h can be regarded 
as adjustable parameters in nuclei. In addition to the 
default parametrization b = h = 0.125, we also use a 
modified one, b = h = 0.06, which gives a smaller ratio 
/ = 1.27 to fit the relative energy between T = 0 and 
T = 1 states in ^°B spectra. 


TABLE I: Adopted parameter sets of effective nuclear inter¬ 
actions. The Bartlett (b) and Heisenberg (h) parameters for 
the Volkov No.2 central interaction and the strength parame¬ 
ter Uls for the GSRS spin-orbit interaction. The Wigner and 
Majorana parameters are fixed to be ui = 0.4 and m = 0.60 
for all sets. The ratio of the to '^E interactions / of the 
central interaction and the Scheerbaum factor Bs of the spin- 
orbit interaction are also shown. 


(A) (B) (A’) (B’) 

b = h 0.125 0.06 0.125 0.06 

/ 1.67 1.27 1.67 1.27 

Uls [MeV] 1300 1300 1000 1000 

Bs [MeV] 103 103 79 79 


III. RESULTS 

We calculate ^°B with the AMD+VAP method. 
AMD wave functions for J'^ = O'*", I’*', 2+, 3"^, 4+ 

states are obtained by VAP. We superpose J^- 
projected states of 10 basis wave functions (five are 
the obtained wave functions and five are the Pp-iSfn- 
projected wave functions) to get energy levels. We 
also apply the AMD-hVAP method to ®Li and ®Be 
and calculate low-lying states, ®Li(l+,2+,3+,0’'') and 
®Be(l/2“, 3/2“, 5/2“, 1/2+, 3/2+, 5/2+). 
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In Fig. [I] we show energy spectra of ®Li and ®Be ob¬ 
tained by the AMD-I-VAP calculation using the inter¬ 
action parameter sets (A) and (B) compared with the 
experimental data. In the ®Li spectra, the level spacing 
between = 1+0, 3+0, and 2+0 states is reproduced 
reasonably. The excitation energy of the 0+1 state is 
overestimated in the result (A) and underestimated in 
the result (B). This means that a value of the ratio / 
in-between / = 1.67 for (A) and / = 1.27 for (B) is rea¬ 
sonable to reproduce the ®Li spectra. It may indicate 
that the effective interaction is slightly weaker in ®Li 
than that in a deuteron. In the ®Be spectra, the excita¬ 
tion energy of the 1/2“ state is reproduced by adjusting 
the spin-orbit strength uis as mentioned previously. Ex¬ 
citation energies of positive-parity states are somewhat 
overestimated, maybe because the present model space 
of AMD wave functions is not sufficient to describe well 
= 1/2+ band states, which are successfully described 
by molecular orbital models [H, . 

We show the calculated energy spectra of ^°B com¬ 
pared with the experimental data in Eig. [51 We also 
show the energy spectra of the NCSM calculation with 
the chiral NN+NNN force. Both results (A) and (B) 
in the present calculation reproduce the ordering of the 
3j^0 and Ij^O states in ^°B. Namely, the 3+0 is the ground 
state and the 1+0 is the first excited state consistently to 
the experimental data and also to the NCSM calculation. 
The relative energy between the 3j^0 and 1+0 states is 
sensitive to the strength of the effective spin-orbit inter¬ 
action. More details of the dependence on the spin-orbit 
interaction and its relation to the NNN force are dis¬ 
cussed later. The Oj^l energy is largely overestimated in 
the result (A) and it is reasonably reproduced in the re¬ 
sult (B) indicating that, in the present model, the smaller 
ratio / ^ 1.27 of the effective ^E and ^E interactions is 
favorable for ^°B than / ^ 1.67 for a deuteron. 

In Table m properties of ®Li, ®Be, and ^°B are listed. 
The present results are compared with the experimental 
data and also theoretical values of the NCSM calculation 
with the chiral NN+NNN force [l^. Properties such as 
radii, moments, and transition strengths are reproduced 
reasonably by the present calculation. 


IV. DISCUSSION 

A. Dependence of energy spectra on spin-orbit 
interaction 

To discuss dependence of the energy spectra on the 
strength of the spin-orbit interaction, we compare the 
energy spectra obtained using the interactions (B) with 
the default strength uis = 1300 MeV and those obtained 
using (B’) with a slightly weak spin-orbit interaction 
uis = 1000 MeV. We show ^°B spectra in Fig. |3l En¬ 
ergies relative to the 3+0 energy of the result (B) are 
plotted. In the result (B’) with a weak spin-orbit inter¬ 
action, the Ij^O energy is lower than the 3j^0 state as 
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FIG. 1: (color online). Energy spectra of ®Li and ®Be ob¬ 
tained by AMD+VAP using the interaction parameter sets 
(A) and (B) compared with the experimental spectra [3^-13^. 
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FIG. 2: (color online). Energy spectra of ^'^B. The theoreti¬ 
cal result of AMD-fVAP using the interaction parameter sets 
(A) and (B), the experimental data [3^. |3^. and the NGSM 
calculation with the chiral NN+NNN force [l^ are shown. 


expected from the pn pair picture that an S = 1 pn pair 
in the 5'-wave is more favored than that in the D-wave 
with no or a weak spin-orbit interaction. As a result, the 
interaction (B’) fails to describe the ordering of low-lying 
energy levels, i.e., the ground state spin, 3+, of ^°B. On 
the other hand, in the result (B), the level inversion of 
the 3)^0 and 1+0 states occurs consistently to the exper¬ 
imental data. The reason for the level inversion is that 
the spin-orbit interaction favors the spin-aligned T = 0 
pn pair and lowers the 3j^0 state, whereas, it gives almost 
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TABLE II: Properties of ®Li, ®Be, and ^°B. Theoretical val¬ 
ues are calculated by AMD-1-VAP using the interactions (A) 
and (B). The experimental proton radii are derived from the 
charge radii in Ref. [^. Other experimental data are taken 
from Refs. [3^-[3^. The values of the NCSM calculation with 
the chiral NN+NNN force from Ref. 0 are also shown. 


Expt. AMD+VAP NCSM 

(A) (B) NN+NNN 


«Li:|A(l+0)| 

31.995 

27.9 

26.4 

32.63 

rp(l^O) [fm] 

2.44(4) 

2.21 

2.21 


Q(l)^O) [e fm^] 

-0.0818(17) 0.09 

0.08 

-0.12(4) 

pil+O) 

0.822 

0.88 

0.88 

0.836 

B{E2-3tO 1^0) 

10.7(8) 

4.3 

4.1 

3.685 

B{E2-,2+0 1+0) 

4.4(23) 

5.5 

5.2 

3.847 

B{Ml-,0tl 1^0) 

15.4(3) 

16.1 

16.4 

15.04(4) 

®Be: 115(3/27)1 

58.164 

53.0 

53.0 


rp(3/27) [fm] 

2.377(12) 

2.42 

2.42 


Q(3/27) [efm^] 

5.288(38) 

5.2 

5.2 


d(3/27) [^iv] 

-1.1778(9) 

-1.24 

-1.24 


i'’B:|£;(3^0)| 

64.751 

58.7 

57.7 

64.78 

rp(3/0) [fm] 

2.28(5) 

2.31 

2.33 

2.197 

Q(3+0) [e fm^] 

8.47(6) 

7.95 

8.2 

6.327 

Ai( 3/'0) [/ijv] 

1.8006 

1.84 

1.85 

1.837 

^(1+0) [^iv] 

0.63(12) 

0.86 

0.84 


R(E2; 1+0 ^ 3)^0) 

4.14(2) 

4.2 

3.6 

3.05(62) 

R(A2;l+0^ 1)^0) 

15.6(17) 

10.2 

10.1 


B(A2;lj0^3^0) 

1.7(2) 

0.9 

1.3 

0.50(50) 

R(A2;2+0 ^ iJO) 

15.2(69) 

2.7 

4.3 


B{E2-,2tO 1^0) 

17.8(18) 

7.9 

7.6 


B(A2;2+0 ^ 3)^0) 

1.2(4) 

1.1 

0.9 


B{E2-,3tO 1^0) 

19.7(17) 

7.7 

8.3 


B{Ml-0+l 1+0) 

7.5(32) 

13.5 

14.7 


B{M1-, itO O/l) 

0.19(2) 

0.0 

0.0 


B{Ml-,2+l 2+0) 

2.5(7) 

3.7 

3.9 


B{Ml-2+l 1+0) 

3.1(8) 

2.8 

2.7 


B{Ml-2+l 1+0) 

0.32(9) 

0.2 

0.4 



no contribution to the energy of the T = 0 pn pair in the 
S'-wave in the 1]*’0 state. 

As for T = 1 states, the 0^1 state somewhat gains 
the spin-orbit interaction energy because the spin-orbit 
interaction favors the T = 1 pn pair in the 0^ state. In 
comparison of the results (B) and (B’), it is found that 
the energy gain for the 0^ state is not as large as that 
for the 3^0 state. 

In general, a spin-orbit interaction in effective two- 
body nuclear interactions used for structure model cal¬ 
culations is an effective spin-orbit interaction in nuclei. 
In principle, it should contain a contribution from the 


TABLE III: Expectation values of harmonic oscillator quanta 
and those of the square spin and angular momentum for ^°B 
calculated with the interaction (A). For harmonic oscillator 
quanta, the minimum value Qmin = 6 for the Ohu configura¬ 
tion is subtracted, and values of AQ — (Q) — Qmin are listed. 


i°B(J’"T) AQ 

(S^) 

(L=) 

CO 

o 

1.0 

2.0 

6.8 

l^O 

1.5 

1.9 

0.1 

O/”! 

0.9 

0.5 

0.5 

1+0 

1.7 

1.9 

5.7 

2+0 

1.4 

2.0 

6.0 

3^0 

1.5 

2.0 

7.0 

2+1 

1.1 

0.5 

6.0 

2+0 

1.5 

2.0 

6.8 

4/-0 

1.1 

2.0 

13.7 
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FIG. 3: (color online). Energy spectra for ^°B calculated with 
AMD-1-VAP using interactions (B) and (B’). The experimen¬ 
tal data are taken from Refs. [3a. |3^ 


NNN force in addition to the original spin-orbit force in 
the bare NN forces. As Kohno pointed out, the NNN 
force contributes attractively to the effective two-body 
spin-orbit interaction in nuclear medium [^ . In the 
G-matrix calculation of nuclear matters using the chiral 
NN and NN + NNN forces in Ref. [2l|, the contribu¬ 
tion of the NNN force is evaluated to be ABs = 20 — 30 
MeV of the Scheerbaum factor. For instance, in a nuclear 
matter with the Fermi momentum kp = 1.35(1.07) fm“^, 
the strength is estimated to be Bs = 84.6(86.5) MeV for 
the chiral NN force and Bs = 116.2(106.7) MeV for the 
chiral NN + NNN force. 

In the present calculation, we use the phenomenologi¬ 
cal effective two-body central and spin-orbit interactions, 
which are adjusted so as to describe the a-a scattering 
and the Is splitting in ®Be. Although the present in¬ 
teractions have no direct link to the bare nuclear forces, 
they should indirectly contain the contribution from the 
NNN force. Then, we reach an idea that a part of the 
two-body spin-orbit interaction in the present effective 
interactions can be interpreted as the contribution from 
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FIG. 4: (color online). Dependence of the energy spectra on 
the strength uis of the spin-orbit interaction for and ®Li 
calculated with AMD-1-VAP. Energy spectra of obtained 
using (B) with the default spin-orbit interaction uis = 1300 
MeV and (B’) with the weaker one uia = 1000 MeV, and those 
of ®Li obtained using (A) with uis = 1300 MeV and (A’) with 
uis = 1000 MeV are shown as well as the experimental en¬ 
ergy spectra. The NCSM calculation using the chiral nuclear 
forces with the NNN force (the chiral NN+NNN force) and 
without the NNN force (the chiral NN force) from Ref. [3 
are also shown. 


the NNN force. With a help of the Kohno’s G-matrix 
analysis, we can roughly estimate the contribution of the 
NNN force in the present parametrization as the change 
Auis ^ 300 MeV of the spin-orbit interaction strength 
which corresponds to the change ABs = 24 MeV of 
the Scheerbaum factor. Therefore, it is expected that 
the result (B’) with the weaker spin-orbit interaction by 
Auis ^ 300 MeV than the default strength can be as¬ 
sociated with the calculation without the NNN force 
contribution in the effective spin-orbit interaction. In 
Fig. m we show energy spectra of calculated with 
the interactions (B) and (B’), and those of the NCSM 
calculations with the chiral NN + NNN and the chiral 
NN forces. In each calculation, the energy of the 3)'’0 
state is set to be zero. As expected, differences in low- 
lying spectra between results (B) and (B’) in the present 
calculation corresponds well to those of the NCSM re¬ 
sults with and without the NNN force, meaning that 
the change Auis ^ 300 MeV of the effective two-body 
spin-orbit interaction gives effects quite similar to the 
contribution of the NNN force on the low-lying spectra 
of ^°B. For instance, the 1)''0 state comes down to the 
lower energy region than the 3)'"0 state in the result (B’) 


because of the reduction Auis ^ 300 MeV consistently 
to the NCSM calculation without the NNN force. The 
excitation energy of the O)*'! state is slightly decreased 
by the reduction Auis ^ 300 MeV, which corresponds to 
the difference of the O)*"! excitation energy between the 
NCSM calculation with the NNN force and that without 
the NNN force. This association of the present results 
(B) and (B’) with the NCSM calculations with and with¬ 
out the NNN force indicates that the part Auis ^ 300 
MeV of the two-body spin-orbit interaction in the present 
phenomenological effective interactions is interpreted as 
the contribution of the NNN force, which is essential 
to the level inversion between the 3)'’0 and 1)*’0 states in 
^°B. We also show ®Li spectra calculated with interac¬ 
tions (A) for the default spin-orbit interaction strength 
and (A’) for the reduced strength, compared with the chi¬ 
ral NN + NNN and NN NCSM calculations. Also for 
®Li, the change in the low-lying spectra by the reduction 
of Auis ^ 300 MeV corresponds well to the difference 
between the NCSM calculations with and without the 
NNN force. 


B. Structure of 

We analyze ^°B wave functions obtained by 
AMD-1-VAP and find that the ground and excited 
states of ^°B are approximately understood by T = 0 
of T = 1 pn pairs around the 2a core. In Table nni we 
show expectation values of the squared intrinsic spin, 
{S'^), and those of the squared orbital angular momen¬ 
tum, (i^). We also show expectation values of the 
harmonic oscillator quanta, (Q), given by the creation 
and annihilation operators Q = a^a of the harmonic 
oscillator for the width parameter v — 0.235 fm“^. Since 
the 2a core gives no contribution to the total intrinsic 
spin, {S'^) reflects mainly intrinsic spin configurations of 
two nucleons around the core. The calculated values of 
for T = 0 states are {S"^) ~ 2 indicating that two 
nucleons form a {ST) = (10) pair, which is the same 
spin-isospin configuration as a deuteron. For T = 1 
states, (5'^) is approximately 0.5 meaning that the T = 1 
pn pair has the dominant {ST) = (01) component with 
a mixing of S' = 1 component. The S = 1 mixing in the 
T = 1 pn pair is nothing but the odd-parity mixing in 
the pair caused by the spin-orbit potential from the core 
as discussed in the previous paper for the pn pair around 
the core in ^®F. {Q) for the 1)''0 state is relatively 
large compared with those for the 3)''0 and O)*"! states 
because the 1)*’0 state has a spatially developed pn pair 
as well as the 2q; clustering and contains higher shell 
components. 

Figure [S] shows the matter density distribution of the 
intrinsic wave functions for the 3)'’0, 1)''0, and O)*'! states. 
The density of the single AMD wave function obtained 
by VAP for each J'^ is shown. In the 3)''0 state, the T = 0 
pn pair exists at the surface of an a cluster, whereas, in 
the 1)''0 state, it spatially develops. In the O)*'! state, the 
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FIG. 5: (color online). Distributions of matter densities 
of '^°B(3)'"), and ^°B(0)*’,T = 1) calculated with 

AMD+VAP using the interaction (A). Densities of intrinsic 
states are integrated with respect to the a axis and plotted 
on the x-y plane (the box size is 10 fm x 10 fm). Here, axes 
of the intrinsic frame are chosen as {x^) > {y'^) > {z^)- 



T = 1 pn pair locates close to an a cluster. As mentioned 
previously, the T = 0 pn pair in the 3)'"0 state and the 
T = 1 pn pair in the O)'’! state are energetically favored 
by the spin-orbit potential from the core. To gain the 
spin-orbit potential, the pn pair remains at the surface 
close to the core in the 3)^0 and O)*"! states. This is 
contrast to the spatially developed T = 0 pn pair in the 
1)*’0 state, in which the spin-orbit interaction gives minor 
contribution. 


V. 2a+pn MODEL ANALYSIS OF pn PAIR 

As discussed previously, the 3)''0 and O)*'! states gain 
the spin-orbit interaction, whereas the 1)'’0 state is not 
affected by the spin-orbit interaction. This result is un¬ 
derstood by effects of the spin-orbit potential to T = 0 
and T = Ipn pairs at the nuclear surface, which were dis¬ 
cussed in the previous paper for based on the ^^0+pn 
model. To reveal the role of the spin-orbit interaction in 
the system, we here apply a 2a+pn model and inves¬ 
tigate effects of the spin-orbit interaction to the pn pair 
at the surface of the 2a core. 

Let us consider a proton and a neutron at the surface 
of the 2a core. Because of the and interactions, 
they form (ST) = (10) and (ST) = (01) pairs. The 
former is the deuteron-like pn pair and the latter corre¬ 
sponds to the dineutron pair. For simplicity, we consider 
two nucleons with parallel intrinsic spins for the T = 0 
pair and untiparallel intrinsic spins for the T = 1 pair 
around the 2a core as shown in Fig. [B] Here we take 
intrinsic spin orientations along the z-axis for the a-a di¬ 
rection. Without the spin-orbit potential from the core, 
it is naively expected that T = 0 and T = 1 pairs move 
in the S'-wave (L = 0) around the 2a core in the lowest 
state to construct J'^T = 1+0 and 0+1 states. Due to 
the stronger ^E interaction than the ^E interaction, the 
1+0 state is expected to be lower than the 0+1 state. In 
the spin-orbit potential from the core, a spin-up nucleon 
at the surface is boosted to have finite momentum and a 
spin-down nucleon is boosted to the opposite direction. 


FIG. 6: Schematic figures for the Sz = 0 pn pair (upper 
figures) and the = 1 pn pair (lower figures) around the 2a 
core in the 2a + pn model. 



(fm) 


FIG. 7: dx dependence of energy calculated with the 2a +pn 
model. Energies of = 0+, 1+ and 3+ projected states 
for ky — 0 are calculated using the interaction (A). The a- 
a distance is fixed to be daa = 3 fm. S'z = 0 is chosen for 
the 0+ state, and = 1 and K = 1{K = 3) are chosen 
for the l+(3+) state. The 3+ energy calculated without the 
spin-orbit interaction is also shown. 


Consequently, for the {ST) = (10) pair, the spin-orbit 
potential boosts two nucleons in the same direction and 
causes the orbital rotation of the pair, and therefore it 
favors the spin-aligned J^T = 3+0 state. For the T = 1 
pair, the spin-orbit potential boosts two nucleons in the 
opposite direction. Due to the opposite boosting by the 
spin-orbit potential, the T = 1 pair is no longer the ideal 
{ST) = (01) pair but it contains the odd-parity mixing, 
i.e., the mixing of the S = 1 component in the dominant 
S = 0 component as discussed in the previous paper. 

To quantitatively discuss contributions of the spin- 
orbit interaction to T = 0 and T = 1 pairs in the 2a+pn 
system, we introduce a, 2a + pn model as follows. The 
2a + pn wave function with anti-parallel spins {Sz = 0) 



















for the T = 1 pn pair is given as 

=A{^o.iRl)<i>c.iR2)MXl)MX2)m) 

1praiX;r) = (/)x{r)Xra, (12) 

where ^a{Rk) is the a cluster wave function written by 
the (Os)'* harmonic oscillator configuration located at Rk, 
and ipra- is the single-particle wave function for a valence 
nucleon assumed to be a localized Gaussian wave packet. 
Here we use labels t = p,n and a =t, i for the isospin 
and intrinsic spin of the nucleon, respectively. We set two 
as with the distance daa parallel to the z-axis as Ri = 
—R 2 = ( 0 , 0 , ^ 0 - 0 / 2 ), and the single-nucleon Gaussian 
wave packets for p f and n | at 


Xi = {dx,iky/‘2v,0), (13) 

X2 = id^,-ikyl2i2,0). (14) 

Here, parameters dx and ky stand for the mean positions 
and momenta of the Gaussian wave packets, 

(</>(Xi, 2 )|f|</)(Xi, 2 )) = (d.,0,0), (15) 

{cl,{x,)\mx^)) = {0,hky,0), (16) 

(</>(X 2 )|p|<^(X 2 )) = (0,-ftfc„0), (17) 


meaning that spin-up and -down nucleons are boosted to 
have finite momenta in the opposite direction (see upper 
panels of Fig. |6l). This parametrization is a kind of ex¬ 
tension of the model for a cluster structures proposed by 
Itagaki et al. in Ref. [i^. Note that, in the ky ^ 0 case, 
the pn pair contains the 5" = 1 component in addition 
to the dominant 5 = 0 component. The 2a + pn wave 
function with parallel spins {Sz = 1) for the T = 0 pn 
pair is written as 

^ttpn = ^{$a(i2l)$a(i?2)V’pt(^l)V'nt(^2)}(4,8) 

with 


Xi = {dx,iky/2n,0), (19) 

X2 = (dx,iky/2n,0), ( 20 ) 

where nucleons in the pn pair are boosted in the same 
direction (see lower panels of Fig. El). 

For simplicity we fix the a-a distance as daa = 3 fm. 
The contribution from the center of mass motion is ex¬ 
actly removed by shifting Gaussian center positions as 
Ri ,2 —^ .R 1.2 — Rg and Xi 2 —^ Xi 2 — Rg with 


Rg = 


4(i2i -|- R2) -t- Xi -|- X2 
10 


( 21 ) 


The J'^ state projected from is given as 


= ( 22 ) 

We calculate energy expectation values of the 2a+pn 
wave functions using the interaction (A) and that with¬ 
out the spin-orbit interaction, and analyze energies of the 
T = 1 and T = 0 pn pairs in the 2a -I- pn system. 


We hrst discuss energies of 2a -\-pn for the ky = 0 case 
with no boosting which corresponds to ideal {ST) — (01) 
and {ST) = (10) pn pairs. Figure [7] shows energies of 
the J'^ = 0+, 1+ and 3+ projected states, 

-Pm® 2 a+pn> and plotted as functions of the 

distance dx of the pair position from the core. Here, 
K = 0, K = 1, and AT = 3 are chosen for = 0+, I’*' 
and 3+ projections, respectively. Note that the = 0+ 
projected wave function is a T = 1 eigen state and 
is a r = 0 eigen state. The J'^ = 0+ and 1+ 
energy curves have minimums in the > 2 fm region 
indicating that the ideal T = 1 and T = 0 pn pairs de¬ 
velop spatially from the core. The spin-orbit interaction 
gives no contribution to the {ST) = (01) pair in the 0+ 
state nor to the {ST) = (10) pair in the I’*' state. In the 
3'*" energy curve obtained without the spin-orbit interac¬ 
tion, the optimum dx at the energy minimum is slightly 
smaller than those for the I’*' and O’*" energy curves be¬ 
cause of the relatively high centrifugal barrier. The 3'*' 
energy obtained with the spin-orbit interaction shows a 
large energy gain in the small dx region. It indicates that 
the spin-aligned T = 0 pair is favored by the spin-orbit 
potential from the core, which keeps the pair close to the 
core. 

Next we analyze the ky ^ 0 case to discuss the contri¬ 
bution of the nucleon momenta. Figure El shows intrinsic 
energies of the 2a -I- pn wave functions for Sz = 0 and 
= 1 without the projection and the J'^ = 0 + pro¬ 
jected energy for Sz = 0. Energies are plotted as func¬ 
tions of the momentum ky. For the J'^ = O'*" projected 
state, we also show (5^), which indicates the 5 = 1 mix¬ 
ing (the odd-parity mixing) in the 5 = 0 component as 
a function of ky. The pn pair position dx is fixed to be 
dx =2 fm. In Figs. ma) and (c) for intrinsic energies, it is 
found that intrinsic states gain the spin-orbit interaction 
in the finite ky region because of the boosting of nucle¬ 
ons in the opposite directions in the 52 = 0 pair and that 
in the same direction in the 5z = 1 pair. In the energy 
curve for the = 0 +-projected state (see Fig. El (b)), a 
further large energy gain of the spin-orbit interaction is 
found in the Hnite ky region. 

In Fig. El we show the O'*' energy with and without the 
spin-orbit interaction plotted on the dx-ky plane. We also 
show the expectation value of the spin-orbit interaction 
of the 0+ projected state. The energy surface obtained 
without the spin-orbit interaction shows the energy min¬ 
imum at dx = 2.2 fm on the ky = 0 line (see Fig. Elb)). 
The contribution of the spin-orbit interaction is attrac¬ 
tive in the Hnite ky, in particular, in the small dx region 
(see Fig. EKc)), in which two nucleons in the 5z = 0 pair 
approximately occupy the single-particle |r2| = |jz| = 3/2 
orbits in the p shell. Gonsequently, the energy minimum 
shifts to the finite ky and slightly smaller dx region in 
the result with the spin-orbit interaction (see Fig. EKa)). 
This result indicates that the {ST) = (01) pn pair in the 
0 "^ state is somewhat broken to contain the odd-parity 
mixing (the 5=1 mixing in the 5 = 0 component) 
by the spin-orbit potential at the surface from the core. 



9 


Moreover, because of the spin-orbit potential, the spatial 
development of the pn pair is suppressed slightly. 

VI. SUMMARY 

We investigated the structures of positive-parity states 
of with AMD-1-VAP using the phenomenological ef¬ 
fective two-body interactions. In the result, we found 
2a + pn structures in We discuss effects of the 
spin-orbit force on the energy spectra and pn correla¬ 
tions in the J^T = 1^0, 3^0, and 0^1 states. The 1^0 
state is not affected by the spin-orbit interaction, whereas 
the 3^0 state gains energy of the spin-orbit interaction 
largely to come down to the ground state, and the 0^1 
state also gains somewhat energy of the spin-orbit inter¬ 
action. We showed that the change Au/g ^ 300 MeV 
of the spin-orbit interaction in the present effective two- 
body interactions gives effects quite similar to the con¬ 
tribution of the NNN force in the NCSM calculation on 
the low-lying spectra of ^°B and ®Li. It indicates that 
the part of the two-body spin-orbit interaction can be 
interpreted as a contribution of the NNN force, which is 
essential to the level ordering of the 3^0 and 1^0 states 
in ^°B. We also applied the 2a + pn model and discuss 
the effects of the spin-orbit interaction on the T = 0 and 
T = 1 pn pairs around the 2a core. In the spin-aligned 
J'^T = 3+0 state, the spin-orbit interaction affects the 
{ST) = (10) pair attractively and suppresses the spa¬ 
tial development of the pair, whereas, in the 1+0 state. 


it gives a minor effect to the {ST) = (10) pair. The 
{ST) = (01) pair in the 0+1 state is somewhat dissoci¬ 
ated to have the odd-parity mixing, i.e., the mixing of 
S' = 1 component by the spin-orbit interaction. 

In the present calculation, we use the phenomenologi¬ 
cal effective two-body central and spin-orbit interactions, 
which are adjusted so as to describe the a-a scattering 
and the Is splitting in ®Be. The present interactions have 
no direct link to the bare nuclear force although the con¬ 
tributions from the NNN force as well as the tensor force 
and also many-body effects in nuclear systems should be 
indirectly contained in the effective interactions. In the 
present paper, we associate a contribution of the NNN 
force with a part of the effective two-body spin-orbit in¬ 
teraction with the help of the G-matrix calculation by 
Kohno. It is a remaining future problem to adopt more 
sophisticated effective interactions derived from bare nu¬ 
clear forces and investigate effects of the NNN force on 
the pn correlations \n Z = N = odd nuclei. It is also 
an important issue to study effects of the NNN force 
on nuclear structures considering the link of the NNN 
force with the effective two-body spin-orbit interactions 
as done for nuclear radii by Nakada et al. [4l|. 
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FIG. 8 : Energies and (S^) calculated with the 2Q;-|-pn model 
using the interaction (A), (a) Energy of the non-projected 
state for the = 0 pn pair, (b) that of the = 0 ^ projected 
state for the Sz = 0 pn pair, and (c) that of the non-projected 
state for the Sz = 1 pn pair. The energies with and without 
the spin-orbit force are plotted as functions of ky. daa = 3 
fm and = 2 fm are chosen, (d) The spin expectation value 
(5'^) of the = 0 ’'" projected state for the Sz = 0 pn pair. 
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FIG. 9: (color online). Energies of the = 0'^ projected 
state for the Sz = 0 pn pair calculated with the 2a + pn 
model using the interaction (A), (a) Energy with the spin- 
orbit interaction, (b) that without the spin-orbit interaction, 
and (c) the expectation value of the spin-orbit interaction. 


















